REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  NO.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions, 
searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments 
regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggesstions  for  reducing  this  burden,  to  Washington 

Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington  VA,  22202-4302. 

Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  oenalty  for  failing  to  comply  with  a  collection  of 
information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY) 

2.  REPORT  TYPE 

3.  DATES  COVERED  (From  -  To) 

New  Reprint 

- 

4.  TITLE  AND  SUBTITLE 

Anomalous  insulator-metal  transition  in  boron  nitride-graphene 
hybrid  atomic  layers 


5a.  CONTRACT  NUMBER 
W911NF-1 1-1-0362 


5b.  GRANT  NUMBER 


6.  AUTHORS 

Li  Song,  Luis  Balicas,  Duncan  Mowbray,  Rodrigo  Capaz,  Kevin  Storr, 
Lijie  Ci,  Deep  Jariwala,  Stefan  Kurth,  Steven  Louie,  Angel  Rubio, 
Pulickel  Ajayan 


5c.  PROGRAM  ELEMENT  NUMBER 

611103 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAMES  AND  ADDRESSES 

William  Marsh  Rice  University 

Office  of  Sponsored  Research 

William  Marsh  Rice  University 

Houston,  TX  77005  - 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND 
ADDRESS(ES) 

U.S.  Army  Research  Office 
P.O.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 
ARO 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

59735 -MS  -MUR.  1 8 


12.  DISTRIBUTION  AVAILIBILITY  STATEMENT 
Approved  for  public  release;  distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 

The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  contrued  as  an  official  Department 
of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other  documentation. 


14.  ABSTRACT 

The  study  of  two-dimensional  (2D)  electronic  systems  is  of  great  fundamental  significance  in  physics.  Atomic 
layers  containing  hybridized  domains  of  graphene  and  hexagonal  boron  nitride  (h-BNC)  constitute  a  new  kind  of 
disordered  2D  electronic  system.  Magnetoelectric  transport  measurements  performed  at  low  temperature  in  vapor 
phase  synthesized  h-BNC  atomic  layers  show  a  clear  and  anomalous  transition  from  an  insulating  to  a  metallic 
behavior  upon  cooling.  The  observed  insulator  to  metal  transition  can  be  modulated  by  electron  and  hole  doping 


15.  SUBJECT  TERMS 

2  DIMENSIONS,  ELECTRONIC- STRUCTURE,  BILAYER  GRAPHENE,  COULOMB  GAP,  BC2N,  RESISTIVITY, 
MANGANITES,  SYSTEMS,  STATE,  FILMS 


16.  SECURITY  CLASSIFICATION  OF: 

a.  REPORT 

UU 

b.  ABSTRACT 

UU 

c.  THIS  PAGE 

UU 

17.  LIMITATION  OF 
ABSTRACT 

UU 


15.  NUMBER  19a.  NAME  OF  RESPONSIBLE  PERSON 

OF  PAGES  Pulickel  Ajayan _ 

19b.  TELEPHONE  NUMBER 
713-348-5904 


Standard  Form  298  (Rev  8/98) 
Prescribed  by  ANSI  Std.  Z39. 1 8 


Report  Title 

Anomalous  insulator-metal  transition  in  boron  nitride-graphene  hybrid  atomic  layers 

ABSTRACT 

The  study  of  two-dimensional  (2D)  electronic  systems  is  of  great  fundamental  significance  in  physics.  Atomic  layers 
containing  hybridized  domains  of  graphene  and  hexagonal  boron  nitride  (h-BNC)  constitute  a  new  kind  of  disordered 
2D  electronic  system.  Magnetoelectric  transport  measurements  performed  at  low  temperature  in  vapor  phase 
synthesized  h-BNC  atomic  layers  show  a  clear  and  anomalous  transition  from  an  insulating  to  a  metallic  behavior 
upon  cooling.  The  observed  insulator  to  metal  transition  can  be  modulated  by  electron  and  hole  doping  and  by  the 
application  of  an  external  magnetic  field.  These  results  supported  by  ab  initio  calculations  suggest  that  this  transition 
in  h-BNC  has  distinctly  different  characteristics  when  compared  to  other  2D  electron  systems  and  is  the  result  of  the 
coexistence  between  two  distinct  mechanisms,  namely,  percolation  through  metallic  graphene  networks  and  hopping 
conduction  between  edge  states  on  randomly  distributed  insulating  h-BN  domains. 


REPORT  DOCUMENTATION  PAGE  (SF298) 
(Continuation  Sheet) 


Continuation  for  Block  13 


ARO  Report  Number  59735.18-MS-MUR 
Anomalous  insulator-metal  transition  in  boron  ni 


Block  13:  Supplementary  Note 

©2012  .  Published  in  Physical  Review  B,  Vol.  Ed.  0  86,  (7)  (2012),  (,  (7).  DoD  Components  reserve  a  royalty-free,  nonexclusive 
and  irrevocable  right  to  reproduce,  publish,  or  otherwise  use  the  work  for  Federal  purposes,  and  to  authroize  others  to  do  so 
(DODGARS  §32.36).  The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  be 
construed  as  an  official  Department  of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other  documentation. 


Approved  for  public  release;  distribution  is  unlimited. 


PHYSICAL  REVIEW  B  86,  075429  (2012) 


Anomalous  insulator-metal  transition  in  boron  nitride-graphene  hybrid  atomic  layers 
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The  study  of  two-dimensional  (2D)  electronic  systems  is  of  great  fundamental  significance  in  physics.  Atomic 
layers  containing  hybridized  domains  of  graphene  and  hexagonal  boron  nitride  (h-BNC)  constitute  a  new  kind  of 
disordered  2D  electronic  system.  Magnetoelectric  transport  measurements  performed  at  low  temperature  in  vapor 
phase  synthesized  h-BNC  atomic  layers  show  a  clear  and  anomalous  transition  from  an  insulating  to  a  metallic 
behavior  upon  cooling.  The  observed  insulator  to  metal  transition  can  be  modulated  by  electron  and  hole  doping 
and  by  the  application  of  an  external  magnetic  held.  These  results  supported  by  ab  initio  calculations  suggest  that 
this  transition  in  h-BNC  has  distinctly  different  characteristics  when  compared  to  other  2D  electron  systems  and 
is  the  result  of  the  coexistence  between  two  distinct  mechanisms,  namely,  percolation  through  metallic  graphene 
networks  and  hopping  conduction  between  edge  states  on  randomly  distributed  insulating  h-BN  domains. 
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I.  INTRODUCTION 

The  interplay  between  Coulomb  interactions  and  random¬ 
ness  has  been  a  long-standing  problem  in  the  physics  of 
2D  electronic  systems.1,2  According  to  the  scaling  theory  of 
localization,  carriers  in  2D  systems  become  localized  due 
to  the  presence  of  randomness  (e.g.,  defects  or  impurities), 
causing  the  electrical  resistance  to  increase  as  the  temperature 
is  lowered  leading  to  an  insulating  ground  state.  Consequently, 
the  nonexistence  of  a  metallic  state  or  the  nonoccurrence  of 
an  insulator  to  metal  transition  (IMT)  in  a  2D  system  was 
predicted.1  Nevertheless,  IMTs  have  been  found  in  a  variety 
of  doped  semiconducting  2DES  (2D-electron  or  2D-hole 
systems),2  and  they  still  represent  an  open  and  intriguing 
problem  for  the  physics  of  2D  systems. 

The  pioneering  work  of  Kravchenko  et  al.  on  the  insulator 
to  metal  transition  in  Si  metal  oxide  semiconductor  field-effect 
transistors  (MOSFETs)3  challenged  the  scaling  theory  of 
localization.1  Since  then,  IMTs  have  been  found  in  several 
semiconducting  2D-electron  or  2D-hole  systems  (2DES).2-4  In 
such  systems,  the  IMT  is  primarily  driven  by  doping;  the  sys¬ 
tem  is  insulating  for  small  carrier  concentrations  and  metallic 
for  large  carrier  concentrations,  with  the  transition  defined  by 
a  critical  carrier  concentration.  The  behavior  of  the  electrical 
resistivity  p  as  a  function  of  temperature  T  characterizes  the 
insulating  <  0)  and  metallic  >  0)  regimes.  Only  for 
a  small  range  of  concentrations  near  criticality  in  the  metallic 
side,  the  system  displays  a  nonmonotonic  behavior  in  p(T). 
Most  recently,  Wang  et  al.  reported  a  dimensionality-driven 
insulator-metal  transition  in  the  Ruddlesden-Popper  series.5 

Pristine  graphene  is  a  2D  allotrope  of  carbon,  and  possesses 
many  unique  electrical  properties.6-8  In  spite  of  its  high  crys¬ 
tallinity,  various  types  of  extrinsic  disorder  such  as  adatoms 


and  charged  impurities,9-11  or  intrinsic  ones  like  surface 
ripples,  topological  defects,  and  the  supporting  substrate  affect 
the  electronic  properties  of  graphene,12-16  which  were  also 
anticipated  by  theoretical  reports.17,18 

Hints  of  IMT  in  graphene  have  been  observed  or  predicted 
recently.19-21  Bolotin  et  al.  have  produced  near-pristine  high- 
mobility  single  layer  graphene  structures,  in  which  ballistic 
conductance  has  been  achieved.19  In  their  devices,  near  the 
neutrality  point,  the  resistance  decreases  with  increasing 
temperature  and  shows  the  opposite  behavior  for  higher 
doping  levels,  increasing  linearly  with  temperature.  Hwang 
and  Das  Sarma20  have  proposed  that  temperature-dependent 
screening  could  explain  these  results.  Finally,  Adam  et  al .21 
found  signatures  of  a  percolation-driven  IMT  in  graphene 
nanoribbbons.  Most  recently,  a  metal/insulator  transition  also 
has  been  reported  in  damaged  graphene  derivatives,  such  as 
oxidized  or  hydrogenated  graphene.22-25 

As  we  shall  see  below,  the  present  results  on  recently 
reported26,27  atomic  layers  containing  hybridized  domains  of 
graphene  and  h-BN  (named  h-BNC)  follow  a  completely 
different  and  novel  phenomenology.  These  systems  have 
already  sparked  a  considerable  amount  of  recent  theoretical 
work.28-31  Specifically,  in  contrast  to  the  above  results  on 
high-purity  or  damaged  graphene  samples  where  the  IMT  is 
externally  imposed  (by  the  gate),  the  IMT  in  h-BNC  is  intrinsic 
and  temperature  dependent,  as  shown  in  the  present  work. 

Here  we  report  electrical  transport  measurements  per¬ 
formed  on  h-BNC  field  effect  transistors  (FET),  and  we  find 
that  the  resistance  of  h-BNC  shows  an  anomalous  insulator 
to  metal  transition  for  all  measured  devices  with  different 
compositions  (from  40  to  94  at.  %  carbon).  With  decreasing 
temperature,  initially  the  resistance  increases  to  a  maximum 
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value  displaying  a  negative  derivative  typical  of  an 
insulating  behavior.  By  further  decreasing  the  temperature,  the 
resistance  drastically  decreases  showing  a  positive  derivative 
typical  of  a  metal. 

This  behavior  is  not  only  robust  under  electron  and  hole 
doping  by  a  gate  potential  and  under  an  external  magnetic 
field,  but  it  also  shows  a  remarkable  single-parameter  scaling 
of  the  R(T)  curves  for  different  doping  and  field  conditions. 
In  particular,  the  doping  dependence  of  R (T)  in  the  metallic 
regime  is  novel,  anomalous,  and,  to  some  extent,  counter¬ 
intuitive:  In  h-BNC,  for  all  at.  %  carbon  compositions,  the 
metallic  regime  is  suppressed,  rather  than  enhanced,  by  doping. 
In  the  insulating  regime,  h-BNC  displays  Efros-Shklovskii 
variable  range  hopping  (ES-VRH)  between  localized 
states.32,33 

The  paper  is  organized  as  follows.  In  Sec.  II  we  describe 
the  methodologies  employed  to  synthesize,  prepare,  char¬ 
acterize,  and  measure  the  h-BNC  devices,  along  with  the 
density  functional  theory  (DFT)  and  nonequilibrium  Green’s 
function  (NEGF)  methods  used  to  model  them.  Resistance 
measurements  for  the  IMT  of  the  h-BNC  devices  are  then 
provided  in  Sec.  Ill,  as  a  function  of  at.  %  carbon  composition, 
backgate  voltage,  and  external  magnetic  field.  These  are 
followed  by  the  results  of  DFT  and  NEGF  calculations  for 
h-BN  nanoribbons  and  nanodomains  embedded  in  graphene. 
In  Sec.  IV  we  discuss  the  observed  single-parameter  scaling 
of  the  IMT  in  h-BNC,  and  combine  our  experimental  and 
theoretical  results  to  develop  a  phenomenological  parallel 
resistor  model  for  the  IMT.  This  is  followed  by  concluding 
remarks. 


II.  METHODOLOGY 

Hybridized  h-BNC  and  pristine  graphene  samples  were 
synthesized  by  thermal  catalytic  chemical  vapor  deposition. 
As  described  in  Ref.  26,  ultrathin  h-BNC  films  were  directly 
grown  on  copper  foil.  After  growth,  PMMA  (poly  methyl 
methacrylate)  was  spin  coated  on  the  Cu  substrates  containing 
grown  h-BNC.  Subsequently,  the  sample  was  immersed  into 
a  diluted  nitric  acid  solution  to  etch  away  the  underlying  Cu 
substrate.  The  h-BNC  films  supported  by  the  PMMA  remained 
floating  on  the  liquid  solution.  The  floating  h-BNC  films  were 
thus  transferred  onto  a  heavily  doped  Si  substrate  having  a 
300  nm  SiC>2  layer,  as  shown  in  Fig.  1(a),  while  PMMA  was 
removed  by  using  acetone  and  isopropanol  solvents. 

To  realize  the  ribbon  geometry,  ultrathin  h-BNC  films 
were  coated  with  the  SI 813  photoresist.  A  UV-mask  aligner 
was  used  to  produce  the  5  to  10  /xm  linear  patterns  on  the 
coated  films  through  optical  lithography.  After  developing, 
the  sample  was  placed  into  an  oxygen  plasma  chamber  to  etch 
away  the  residual  parts,  leaving  only  the  ribbons,  as  shown 
in  Fig.  1(b).  After  removing  the  photoresist  with  acetone,  the 
h-BNC  and  pristine  graphene  ribbons  with  length  10-30  /xm 
were  contacted  with  Ti/Au  (3/30  nm)  metal  electrodes.  The 
Ti/Au  metallization  was  used  to  create  Ohmic  contacts  on  the 
h-BNC  ribbon  by  optical  lithography  and  £-beam  evaporation. 
Examples  of  the  resulting  h-BNC  devices  are  seen  in  the  optical 
and  scanning  electron  microscopy  images  of  Figs.  1(c)  and  2, 
respectively. 


FIG.  1 .  (Color  online)  Optical  microscopy  images  of  a  transferred 
h-BNC  film  (a)  and  ribbon  (b)  on  a  silicon  substrate.  The  patterned 
ribbon  was  fabricated  by  using  optical  lithography  and  oxygen  plasma 
etching,  (c)  Optical  image  of  a  typical  h-BNC  ribbon  electrically 
contacted  with  Ti/Au  pads  of  thickness  ranging  from  3  to  30  nm. 


Figure  3  shows  a  schematic  depicting  the  geometry  of  the  h- 
BNC  devices.  The  drain  to  source  voltage  was  applied  through 
the  top  two  electrodes  while  the  backgate  voltage  was  loaded 
from  the  bottom  Si++  electrode. 

Figure  4(a)  shows  a  typical  high-resolution  transmission 
electron  microscopy  (HRTEM)  image  of  the  edge  of  our  tested 
films.  It  can  be  clearly  seen  that  our  h-BNC  films  mostly  consist 
of  three  layers.  The  fast  Fourier  transform  (FFT)  pattern  in 
Fig.  4(b)  taken  from  the  plane  of  the  films  shows  three  sets  of 
hexagonal  spots  with  a  rotational  angle  of  1 1°,  which  indicates 
that  our  h-BNC  is  a  three-layer  hexagonal  packing  region. 
The  electron  energy  loss  spectroscopy  (EELS)  performed  on 
the  h-BNC  film  to  determine  their  chemical  composition  and 
structure  is  shown  in  Fig.  4(c).  This  proves  that  the  films 
consist  of  B,  C,  and  N.  These  three  elements  are  sp2  hybridized 
and  uniformly  distributed  throughout  the  film.  More  details  of 
the  structure  characterization  of  h-BNC  can  be  found  in  our 
previous  work.26 


FIG.  2.  Scanning  electron  microscope  image  of  one  of  our 
experimental  devices. 
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FIG.  3.  (Color  online)  Schematic  diagram  of  the  hybridized  h-BN 
graphene  device  geometry. 


We  note  that  all  the  experimental  data  presented  here  are 
based  on  samples  composed  of  two  to  three  layers  instead  of 
a  single  layer  of  h-BNC.  All  the  subsequent  low-temperature 
transport  measurements  were  performed  by  using  a  variable 
temperature  (0.5-300  K)  3 He  cryostat  under  magnetic  fields 
up  to  9  T.  A  variety  of  ac  and  dc  setups  were  employed,  with 
the  external  magnetic  field  ( B )  applied  normal  to  the  h-BNC 
plane. 


FIG.  4.  (a)  High-resolution  transmission  electron  microscopy 
(HRTEM)  image  taken  on  the  edge  of  the  sample  indicates  most 
of  our  tested  films  are  composed  of  three  layers,  (b)  The  fast  Fourier 
transform  reveals  a  three-layer  packing  region  with  rotational  angle 
of  1 1°  in  our  sample,  (c)  The  core  electron  energy  loss  spectra  (EELS) 
clearly  shows  B,  C,  and  N  are  uniformly  distributed  throughout  the 
film. 


For  the  low-temperature  transport  measurements,  h-BNC 
samples  were  deposited  on  poly  silicon  hetero  structures,  typ¬ 
ically  covered  with  a  300  nm  thick  insulating  SiC>2  layer 
which  prevents  any  leakage  of  carriers  from  the  sample 
to  the  substrate.  Samples  were  measured  using  the  source- 
measurement  unit  Keithley  2602A  (which  has  a  sensor  input 
impedance  >10  GQ).  For  samples  with  lower  resistances, 
we  checked  the  obtained  results  by  using  either  a  Stan¬ 
ford  Research  System  830  Lock-In  (input  impedance  ~10 
MQ)  in  combination  with  a  Keithley  6221  current  source 
(output  resistance:  >100  TQ  in  the  2  to  20  nA  range) 
or  a  Lakeshore  370  resistance  bridge  (input  impedance 
approaching  5  TQ)  with  excitation  currents  ranging  from  0.01 
to  1  nA. 

We  also  checked  that  the  obtained  results  were  current 
independent.  Samples  were  mounted  onto  male  chip  carriers, 
themselves  plugged  into  female  chip  carriers  adapted  to 
rotators  and  wired  with  polyimide  insulated  phosphor  bronze 
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FIG.  5.  (Color  online)  Low-temperature  transport  of  h-BNC  and 
pristine  graphene  FETs.  Temperature  dependence  of  the  resistance 
of  a  pristine  graphene  (a)  and  a  hybridized  BNC  with  90  at.  % 
carbon  (b)  under  zero-magnetic  field  in  a  cooling  run.  The  curve 
of  h-BNC  exhibits  a  peak  at  rpeak  =  12  K,  signaling  the  insulator- 
metal  transition.  A  negative  derivative  of  the  resistance  (ff  <  0) 
corresponds  to  an  insulating  phase,  and  a  positive  sign  (§  >  0)  is 
empirically  associated  with  the  metallic  phase.  It  should  be  clarified 
that  graphene  was  measured  in  a  four  terminal  configuration  while 
h-BNC  was  measured  in  a  two  terminal  one.  The  insert  images  of 
(a)  and  (b)  are  atomic  models  for  pure  graphene  and  hybridized  h-BN 
graphene,  respectively. 
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FIG.  6.  (Color  online)  Temperature  dependence  of  the  resistance  for  h-BNC  films  with  (a)  90  at.  %  carbon  under  zero  gate  voltage  and 
under  various  magnetic  fields  measured  either  in  warming  up  or  cooling  down  sweeps.  Insulator  to  metal  transition  observed  in  different  h-BNC 
devices  with  (b)  65  at.  %,  (c)  94  at.  %,  and  (d)  56  at.  %  carbon. 


twisted  pairs  where  each  wire  was  soldered  to  a  pin  whose 
typical  separation  is  of  nearly  1  mm  and  which  were  further 
insulated  with  heat  shrinking  polyolefin  tubing.  3 He  was 
used  as  an  exchange  gas.  We  have  measured  more  than  five 
h-BNC  devices  with  different  compositions  (the  percentage  of 
carbon  was  defined  by  XPS)  in  order  to  confirm  this  unique 
insulator-metal  transition,  and  each  device  has  also  been  tested 
several  times  to  avoid  any  artifacts. 

All  DFT  calculations  have  been  performed  using  the  grid- 
based  projector  augmented  wave  function  code  GPAW, 34,35 
employing  a  double-^  -polarized  locally  centered  atomic  or¬ 
bital  (LCAO)  basis  set,  and  the  Perdew-Burke-Ernzerhof 
(PBE)  exchange-correlation  functional.36  The  occupation  of 
the  Kohn-Sham  orbitals  was  calculated  at  an  electronic  tem¬ 
perature  of  kBT  ~  0.1  eV,  with  all  energies  then  extrapolated 
to  T  =  0  K.  The  product  of  the  supercell  dimensions  and 
the  number  of  k  points  used  in  the  calculations  is  >25  A 
in  all  repeated  directions.  Spin  polarization  has  been  taken 
into  account  in  all  calculations.  Structural  relaxation  has  been 
performed  until  a  maximum  force  below  0.05  eV/A  was 


reached.  Transport  calculations  have  been  performed  using 
the  NEGF  method  with  the  electronic  Hamiltonian  obtained 
from  the  DFT  calculation. 

III.  RESULTS 

Figure  5(a)  shows  the  temperature  dependence  of  the 
resistance  for  a  pristine  graphene  sample  in  zero-magnetic 
field.  To  provide  the  resistivity  instead  of  the  raw  resistance 
traces,  one  would  need  to  know  the  precise  number  of  BNC 
layers  in  the  device,  which  in  our  case  could  not  be  determined. 
As  is  clearly  seen,  the  resistance  grows  monotonically  as  the 
temperature  decreases,  which  indicates  that  pure  graphene  acts 
as  a  semiconductor  with  zero  band  gap  in  agreement  with 
previous  reports.6-8 

Figure  5(b)  shows  a  typical  temperature  dependence  of  the 
resistance  at  zero  field  for  an  h-BNC  sample  with  90  at.  % 
carbon.  In  contrast  with  pristine  graphene,  h-BNC  shows  a 
clear  insulator  to  metal  transition  with  a  peak  at  around  12  K. 
Here  a  negative  derivative  of  the  resistance  as  a  function  of 
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temperature  (^)  corresponds  to  an  insulating  behavior  and  a 
positive  sign  is  empirically  associated  with  a  metallic  behavior, 
as  marked  in  Fig.  5(b).  Experimentally,  on  the  insulating 
side  of  the  transition,  the  two  terminal  resistance  of  h-BNC 
monotonically  increases  to  reach  a  maximum  at  ~4837  k£2 
as  the  temperature  decreases  from  room  temperature  to  12  K. 
On  the  metallic  side  of  the  transition,  the  resistance  drops 
drastically  by  more  than  one  order  of  magnitude  to  ~  193  k£2 
as  the  temperature  decreases  from  12  K  to  the  lowest  achieved 
value  of  500  mK. 

We  have  measured  our  samples  in  constant  electrical 
current  mode,  with  currents  ranging  from  0.01  to  1  nA,  not 
observing  any  significant  effect  of  the  magnitude  of  the  current 
across  the  transition,  meaning  no  significant  differences  in  the 
respective  values  of  the  resistivity  in  any  range  of  temperatures. 
We  remark  that  all  the  measurements  are  within  the  Ohmic 
regime  in  the  transition  region. 

The  same  R  versus  T  curves  for  h-BNC  devices  with  other 
at.  %  carbon  compositions  are  also  provided  in  Figs.  6  and  7. 
Figure  6(a)  shows  a  typical  resistivity  measurement  under  zero 
gate  electric  field  and  under  magnetic  fields  up  to  9  T  for  the 
h-BNC  sample  with  90  at.  %  carbon.  It  indicates  from  the 
warming  and  cooling  down  temperature  sweeps  that  a  very 
small  hysteresis  is  observed  between  temperature  up  and  down 
sweeps. 

Figures  6(b)-6(d),  and  7  show  several  typical  traces  of  the 
resistance  as  a  function  of  temperature  for  h-BNC  films  with 
different  carbon  aspect  ratios,  including  65,  94,  56,  and  40 
at.  %  carbon  percentages,  respectively.  We  have  observed  the 
transition  in  all  of  our  samples  except  for  the  pure  material 
(C  or  BN).  These  devices  were  measured  at  two  independent 
laboratories  with  different  experimental  setups,  and  in  both 
circumstances  the  same  behavior,  namely  an  insulator  to  metal 
transition,  was  observed  in  our  samples.  Notice  that  even 
samples  displaying  much  lower  resistances  (nearly  two  orders 
of  magnitude  smaller),  such  as  the  device  fabricated  from 
BNC  containing  94  at.  %  carbon  percentage,  also  display 
the  same  transition  (Fig.  7).  This,  by  itself,  is  enough  to 


FIG.  7.  (Color  online)  Temperature  dependence  of  the  resistance 
of  a  h-BNC  film  under  gate  voltages  of  Vgate  =  +40,  0,  —40  V,  and 
under  fixed  magnetic  field  of  0  T  for  a  40  at.  %  carbon  device. 


Temperature  (K) 


Temperature  (K) 

FIG.  8.  (Color  online)  Temperature  dependence  of  the  resistance 
of  a  h-BNC  film  under  gate  voltages  of  Vgate  =  +40,  +28,  0,  —40  V, 
and  under  fixed  magnetic  field  of  (a)  6  T  and  (b)  9  T  for  a  90  at.  % 
carbon  device. 


discard  any  instrumental  artifact,  such  as  leakage  currents  or 
the  proximity  of  the  sample’s  resistance  to  the  input  impedance 
of  the  measuring  device. 

Figures  7  and  8  show  the  temperature  dependence  of  the 
resistance  of  a  h-BNC  film  for  various  applied  backgate 
voltages  under  fixed  magnetic  fields  of  0  T  for  a  40  at.  % 
carbon  device,  and  6  and  9  T  for  a  90  at.  %  carbon  device, 
respectively.  It  is  clearly  seen  that  the  peak  position  defining 
the  IMT  moves  to  lower  temperatures  with  increasing  backgate 
voltage  for  either  electrons  (within  the  range  of  —40  to  28  V)  or 
holes  (within  the  range  of  +28  to  +40  V).  This  peak  becomes 
sharper  as  it  moves  to  lower  temperatures. 

The  changes  in  resistance  as  a  function  of  backgate  voltage 
were  measured  at  a  temperature  of  37.5  K  on  the  h-BNC  FET 
devices.  As  is  clearly  shown  in  Fig.  9(b),  the  resistance  curves 
exhibit  an  ambipolar  semiconducting  behavior  for  h-BNC, 
which  is  similarly  observed  in  pristine  graphene.  Considering 
the  typical  p- type  semiconducting  behavior  of  B-C-N  mixed 
atomic  structures  as  expected  from  both  theoretical  predictions 
and  previously  reported  data,37-39  we  propose  that  the  behavior 
of  our  h-BNC  films  can  be  understood  as  a  percolating 
graphene  network  with  embedded  BN  domains.26,27  Note  that 
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Drain  Voltage  (V) 

FIG.  9.  (Color  online)  Electric  field  dependence  of  the  insulator-metal  transition  in  h-BNC  at  90  at.  %  carbon,  (a)  Temperature 
dependence  of  the  h-BNCs  resistance  at  gate  voltage  Vgate  =  +40,  +28,  0,  —28,  —40  V  in  a  fixed  magnetic  field  of  3  T.  (b)  The  variation 
of  resistance  as  a  function  of  gate  voltage  measured  at  a  temperature  of  37.5  K.  The  vertical  line  marks  the  neutrality  point  (Dirac 
point)  of  h-BNC  FET.  (c)  The  corresponding  changes  of  peak  position  and  FWHM  in  abnormal  IMT  peak  as  a  function  of  gate  voltage 
in  (a),  (d)  Room  temperature  I-V  curves  for  h-BNC  under  different  bias  voltages. 


the  maximum  resistance  value  (corresponding  to  the  Dirac 
point)  for  h-BNC  is  shifted  to  28  V  as  a  result  of  the  chemical 
doping  and  boundary  scattering  between  h-BN  and  graphene 
domains.  The  calculated  electron  mobility  of  our  h-BNC  de¬ 
vices  is  12  cm2  V-1  s-1  with  a  hole-mobility  of  8  cm2  V-1  s_1, 
which  is  much  smaller  than  the  mobilities  reported  for 
graphene.6-8  This  indicates  that,  although  the  sp 2  network 
is  preserved  in  h-BNC,26,2  the  random  distribution  of  h-BN 
clusters  and  associated  domain  boundaries  are  strong  sources 
of  electron  scattering,  typical  of  highly  disordered  systems. 

We  also  analyzed  the  IMT  in  h-BNC  as  a  function  of 
doping  by  measuring  the  R(T )  curves  while  modulating  the 
gate  voltage  within  the  range  of  —40  to  40  V,  as  shown 
in  Fig.  9.  Figure  9(a)  shows  the  temperature  dependence 
of  the  resistance  of  h-BNC  at  various  doping  levels  in  a 
fixed  magnetic  field  of  3  T.  We  observe  in  Fig.  9(b)  that  the 
conductance  of  h-BNC  becomes  a  -dominated  conductance 
within  the  gate  voltage  range  (—40  <  Vgate  <  28  V),  while 
it  shows  an  ft -dominated  conductance  when  the  gate  voltage 
is  larger  than  28  V.  Here  we  found  that  the  IMT  transition 
temperature  for  h-BNC  with  -dominated  conductance  down¬ 


shifts  from  10.5  to  3  K  when  the  gate  voltage  is  varied  from 
28  to  —40  V,  while  the  IMT  transition  temperature  of  h-BNC 
with  ft -dominated  conductance  drops  down  to  9.5  K  when  the 
gate  voltage  is  increased  from  28  to  40  V.  Meanwhile,  the 
full  width  at  half  maximum  (FWHM)  of  the  peak  observed 
at  the  IMT  decreases  as  the  peak  position  is  shifted  to  lower 
temperatures. 

Figure  9(c)  shows  the  changes  in  position  and  FWHM  of 
the  IMT  peak  as  a  function  of  the  gate  voltage  when  measured 
in  a  magnetic  field  of  3  T.  It  indicates  that  the  IMT  peak 
down-shifts  and  its  FWHM  decreases  with  increasing  doping 
for  both  electrons  and  holes. 

Figure  9(d)  reproduces  a  series  of  I-V  characteristics  for 
the  h-BNC  sample  with  90  at.  %  carbon  at  room  temperature 
under  different  bias  voltages.  The  linear  behavior  seen  here 
indicates  that  the  contact  resistance  between  h-BNC  samples 
and  the  metallic  electrodes  in  our  experiment  is  not  significant. 
Moreover,  by  comparing  the  data  for  the  two-probe  and  four- 
probe  configurations  from  our  previous  work,26  we  found  no 
significant  differences  between  the  two  types  of  conductivity 
data  obtained  for  our  devices. 
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FIG.  10.  (Color  online)  Temperature  dependence  of  the  resistance 
of  a  h-BNC  film  under  magnetic  fields  of  B  =  0,  3,  6,  and  9  T, 
measured  for  backgate  voltage  values  of  (a)  —40  V  and  (b)  +40  V 
for  a  90  at.  %  carbon  device. 


Figure  10  shows  the  IMT  behavior  of  h-BNC  as  a  function 
of  magnetic  field  for  fixed  values  of  backgate  voltages  ranging 
from  —40  to  +40  V,  respectively.  It  indicates  that  the  IMT 
transition  peak  shifts  to  lower  temperatures  and  becomes 
narrower  when  an  external  magnetic  field  is  applied.  The 
detailed  dependence  of  the  IMT  observed  in  h-BNC  on  either 
magnetic  or  electric  field  is  shown  in  Table  I.  The  lowest 
observed  transition  temperature  is  2.95  K  under  an  applied 
field  of  9  T  and  with  an  applied  backgate  voltage  of  —40  V  on 
the  h-BNC  film. 

We  measured  the  anomalous  IMT  in  h-BNC  as  a  function 
of  perpendicular  magnetic  fields  up  to  9  T.  Figure  1 1(a)  shows 
the  temperature  dependence  of  the  resistance  of  h-BNC  under 
several  magnetic  field  values  for  zero  gate  voltage.  As  seen,  the 
IMT  transition  temperature  shifts  to  lower  temperatures  when 
an  external  magnetic  field  is  applied.  The  higher  the  magnetic 
field  the  larger  the  observed  down-shift.  Notice  also  that  the 
peak  at  the  IMT  becomes  narrower  and  sharper  as  the  peak 
position  is  down- shifted  to  lower  temperatures.  The  FWHM 
of  the  IMT  peak  decreases  with  increasing  magnetic  field. 


TABLE  I.  Comparison  of  insulator-metal  transition  peak  position 
Fpeak  and  full  width  at  half  maximum  (FWHM)  in  K  of  h-BNC  with  90 
at.  %  carbon  at  various  magnetic  field  strengths  B  in  T  and  backgate 
voltages  Vgate  in  V. 


Vgate  (V) 

B(  T) 

0 

3 

6 

9 

40 

Tpeak  (K) 

10.90 

9.43 

9.27 

8.15 

FWHM  (K) 

12.50 

11.66 

10.32 

9.28 

28 

Tpeak  (K) 

11.50 

10.75 

9.93 

9.30 

FWHM  (K) 

14.30 

12.80 

11.57 

10.51 

0 

Fpeak  (K) 

10.38 

9.18 

8.72 

7.79 

FWHM  (K) 

13.10 

12.30 

10.60 

9.14 

-40 

Tpeak  (K) 

3.70 

3.06 

2.98 

2.95 

FWHM  (K) 

4.87 

4.17 

3.55 

3.28 

Figure  11(b)  indicates  the  changes  in  position  and  FWHM 
of  the  IMT  peak  as  a  function  of  magnetic  field  measured 
under  a  gate  voltage  of  28  V.  As  clearly  seen  by  the  position 
of  the  peak,  the  FWHM  drops  monotonically  with  increasing 
magnetic  field.  The  dependence  of  IMT  on  magnetic  field 
measured  for  various  gate  voltages  is  listed  in  Table  I.  These 
observations  on  the  IMT  and  its  dependence  on  gate  voltage 
and  magnetic  field  are  reproducible  and  qualitatively  similar 
for  different  h-BNC  samples  having  different  carbon  content. 
We  remark  that  in  contrast  to  2DEG  (MOSFETS)4  where  the  B 
field  completely  suppresses  the  IMT  transition,  in  the  present 
case  for  h-BNC  the  IMT  transition  shifts  to  lower  temperatures 
but  remains  robust. 

To  better  understand  the  underlying  mechanisms  behind 
the  IMT  in  h-BNC  systems,  we  have  performed  ab  initio 
calculations  at  the  DFT  level  for  several  model  h-BNC  type 
systems.  These  include  h-BN  nanoribbbons  of  various  widths 
embedded  in  graphene,  depicted  schematically  in  Fig.  12,  and 
h-BN  nanodomains  embedded  in  graphene,  shown  in  Fig.  13. 

Figures  12  and  13(a)  depict  Kohn-Sham  states  localized  on 
the  edges  of  h-BN  nanoribbons  and  nanodomains  embedded 
in  graphene,  respectively.  For  both  systems,  these  states  are 
near  the  Fermi  level  (at  0  and  0.2  eV,  respectively),  suggesting 
that  they  may  be  the  states  through  which  ES-VRH  occurs. 

Transport  calculations  have  been  performed  using  the 
NEGF  method  with  the  electronic  Hamiltonian  obtained  from 
the  DFT  calculation.  The  calculated  transmission  through  a 
N  terminated  zigzag  edged  h-BN  nanodomain  embedded  in 
graphene  is  shown  in  Fig.  14,  with  a  structural  schematic  as  an 
inset.  The  zero-bias  transmission  at  the  Fermi  level  is  ~0.2, 
giving  a  resistance  of  ~60  kQ.  This  agrees  qualitatively  with 
the  low-temperature  limit  of  resistance  in  the  experimental 
system  of  200  kQ.  It  should  be  noted  that  the  peak  in 
transmission  at  ~0.2  eV  corresponds  to  the  position  of 
the  state  localized  on  the  N  terminated  zigzag  edge  of  the 
h-BN  nanodomain,  shown  in  Fig.  13(a).  Calculations  of  the 
transmission  across  h-BN  nanoribbons  embedded  in  graphene, 
shown  schematically  in  Fig.  12,  indicate  these  structures  act 
as  insulators,  with  resistances  in  the  TQ  range,  even  for  the 
smallest  nanoribbon  widths  considered. 

Investigations  have  also  been  performed  into  how  doping 
via  an  external  gate  voltage  influences  the  states  in  h-BNC 
systems.  This  has  been  accomplished  by  performing  charged 
calculations  for  both  zigzag  terminated  h-BN  nanoribbons  and 
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FIG.  1 1 .  (Color  online)  Magnetic  field  dependence  of  the  insulator-metal  transition  in  h-BNC  at  90  at.  %  carbon,  (a)  Temperature  dependence 
of  the  h-BNCs  resistance  at  magnetic  field  B  =  0,  3,  6,  and  9  T  measured  in  zero  gate  voltage.  The  abnormal  IMT  point  down-shifts  to  lower 
temperature  with  the  increasing  magnetic  field.  The  dashed  lines  are  guides  to  the  eye.  (b)  Changes  of  peak  position  and  FWHM  in  the  abnormal 
IMT  peak  as  a  function  of  magnetic  field  measured  in  a  gate  voltage  of  28  V. 


N-terminated  zigzag  edged  h-BN  nanodomains  embedded  in 
graphene.  In  this  case,  a  single  valence  electron  has  been 
removed  from  the  supercell,  and  replaced  by  a  uniform 
background  electric  potential  to  ensure  no  dipole-dipole  inter¬ 
action  occurs  between  neighboring  supercells  in  the  periodic 
calculation.  Such  calculations  effectively  mimic  the  effect  of 
an  external  gate  voltage  (or  doping)  on  the  system. 

Figure  15  depicts  structural  schematics  and  the  change  in 
charge  density  A p  upon  removal  of  one  electron  for  h-BN 
zigzag  terminated  nanoribbons  and  N  terminated  zigzag  edged 
h-BN  nanodomains  embedded  in  graphene.  While  charge  is  re¬ 
moved  from  the  extended  graphenelike  states  in  both  systems, 
it  is  also  transferred  from  the  localized  edge  states  on  the  h-BN 
nanoribbons  and  nanodomains  to  more  tightly  bound  states 
also  located  on  these  edges.  More  importantly,  the  energy  of  the 
localized  edge  states  is  nearly  unchanged  relative  to  the  Fermi 


i  w  w  w  w  w  | ' 

wjf  »  W 

w 

V 

V 

©  ©  ©  ©  ©  © 

>  ©  ©  ©  ©  Q© 

-  ©  GO 

o 

© 

©  ©  ©  ©  ©  ©  Km 

t  ©  ©  ©  ©  Q©  1 

go 

o 

© 

©  ©  ©  ©  ©  © 

>  ©  ©  ©  ©  Q©  1 

©  GO 

o 

© 

©  ©  ©  ©  ©  © 

t  ©  ©  ©  ©  Q©  p 

©  GO 

o 

© 

©  ©  ©  ©  ©  © 

»  ©  ©  ©  ©  Q©  p 

^  GO 

o 

© 

©  ©  ©  ©  ©  ©  Qy 

>  ©  ©  ©  ©  Q  ©  p 

GO 

o 

© 

©  ©  ©  ©  ©  ©  Qq 

*  ©  ©  ©  ©  Q  ©  p 

GO 

o 

© 

©  ©  ©  ©  ©  © 

>  ©  ©  ©  ©  Q©  p 

GO 

o 

© 

©  ©  ©  ©  ©  ©  K4© 

t  ©  ©  ©  ©  Q©  p 

^  •• 

e 

© 

©  ©  ©  ©  ©  ©  U  G 

>  ©  ©  ©  ©  Q©  p 

Lo  go 

o 

© 

©  ©  ©  ©  ©  © 

©*©  G 

G  o 

>  ©©  ©©  Q©  p 

GO 

o 

© 

©  ©  ©  ©  ©  ©  ft^G 

*  ©  ©  ©  ©  Q©  p 

^  GO 

o 

© 

©  ©  ©  ©  ©  ©  U  G 

>  ©  ©  ©  ©  Q©  p 

Lo  go 

o 

© 

A 

A 

A 

V#  w  w  \ 

©  ©  ©  © 

©  G  G  4 

GO  GO 
G  G  G  4 

G  G  G  G 

G  G  G  4 

G  G  G  G 

G  G  G  4 

G  G  G  G 

G  G  G  4 

G  G  G  G 

G  G  G  4 

G  G  G  G 

G  G  G  4 

G  G  G  G 

G  G  G  4 

G  G  G  G 

G  G  G  4 

G  G  G  G 

G  G  G  4 

G  G  G  G 

G  G  G  4 

G  G  G  G 

G  G  G  4 

G  G  G  G 

G  G  G  4 

fi  lA  tj  ti 


level.  This  suggests  that  for  the  level  of  doping  considered  in 
the  experiments  discussed,  these  localized  edge  states  are  kept 
near  the  Fermi  level,  and  available  for  VRH  to  occur. 


IV.  DISCUSSION  AND  CONCLUSIONS 

Figure  9(a)  shows  that  the  peak  and  minimum  values  for 
the  resistance  as  a  function  of  temperature  are  remarkably 
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FIG.  12.  (Color  online)  Schematic  of  a  h-BN  nanoribbon  em¬ 
bedded  in  graphene.  Two  degenerate  states  localized  at  the  N  and 
B  terminated  zigzag  edges  of  the  h-BN  nanoribbon  are  depicted  by 
isosurfaces  of  ±0.1  e/A3/2.  B,  N,  and  C  atoms  are  depicted  by  red, 
blue,  and  gray  balls,  respectively. 


FIG.  13.  (Color  online)  Schematics  of  a  N  terminated  zigzag 
edged  nanodomain  in  h-BNC  with  wave  functions  depicted  as 
isosurfaces  of  ±0.025  A-3/2  for  (a)  localized  edge  states  at  the 
nanodomain  boundaries  and  (b)  extended  graphenelike  states.  B,  N, 
and  C  atoms  are  colored  red,  blue,  and  gray,  respectively. 
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FIG.  14.  (Color  online)  Transmission  through  a  N  terminated 
zigzag  edged  h-BN  nanodomain  embedded  in  graphene  (shown  as 
inset)  as  a  function  of  energy  s  relative  to  the  Fermi  level  sF  in 
eV.  Note  the  peak  in  the  transmission  at  ~0.2  eV  corresponding  to 
the  state  localized  on  the  edge  of  the  h-BN  nanodomain,  shown  in 
Fig.  13(a). 


independent  of  the  doping  levels  for  the  90  at.  %  carbon 
sample.  This  suggests  that  doping  and  temperature  effects 
must  be  described  in  a  combined  way  and  that  the  R(T)  curves 
should  display  a  single-parameter  scaling.  In  other  words,  all 
resistance  curves  should  follow  a  universal  function  of  the 
rescaled  temperature  T /  rpeak-  Single-parameter  scaling  is  a 
signature  feature  of  IMT  in  semiconducting  2DES. 
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FIG.  15.  (Color  online)  Schematic  of  a  (a)  h-BN  zigzag  ter¬ 
minated  nanoribbon  and  a  (b)  N  terminated  zigzag  edged  h-BN 
nanodomain  embedded  in  graphene.  The  change  in  charge  density 
A p  upon  removal  of  one  electron  is  shown  by  isosurfaces  of 
±0.00625  e/A3.  Red  regions  depict  positive  charge,  while  blue 
regions  depict  negative  charge.  B,  N,  and  C  atoms  are  depicted  by 
red,  blue,  and  gray  balls,  respectively. 


(a)  Ftm(T) 


FIG.  16.  (Color  online)  Theoretical  model  for  the  insulator-metal 
transition  in  h-BNC:  (a)  Phenomenological  parallel-resistor  model 
for  the  IMT.  (b)  Single-parameter  scaling  leading  to  the  collapse  of 
all  five  R(T)  curves  shown  in  Fig.  9(a),  for  different  doping  levels, 
into  a  single  universal  curve.  The  parallel-resistor  model  (black  line) 
is  compared  to  the  experimental  resistance  for  the  five  different  gate 
voltages. 

Figure  16  shows  that  indeed  all  R(T )  curves  for  gate 
voltages  —40,  —28,  0,  ±28,  and  ±40  V  collapse  into  a 
single  universal  curve  when  plotted  as  a  function  of  the 
rescaled  temperature.  The  collapse  is  nearly  perfect,  which 
is  an  indication  of  the  similarity  between  the  IMT  in  h-BCN 
and  in  2DES.  Remarkably,  the  scaling  works  equally  well  in 
both  the  metallic  and  the  insulating  regions. 

On  the  other  hand,  Fig.  7  shows  that  the  peak  resistance 
for  the  40  at.  %  sample  is  doping  dependent,  which  implies 
the  absence  of  single-parameter  scaling,  even  though  the  IMT 
is  very  clear  in  this  sample.  This  may  be  an  indication  that 
IMT  for  samples  with  low  carbon  content  may  follow  a 
slightly  different  physics  than  the  IMT  for  samples  with  high 
carbon  content.  Indeed,  it  seems  reasonable  that  a  percolating 
graphene  network  with  embedded  “insulating”  BN  clusters 
should  have  different  electrical  transport  properties  than  a 
percolating  BN  network  embedded  with  “metallic”  graphene 
clusters.  Further  experimental  and  theoretical  work  is  needed 
to  address  this  issue. 

In  any  case,  there  are  important  and  surprising  differences 
between  the  IMT  in  h-BNC  and  2DES.  First,  as  we  mentioned, 
the  IMT  in  h-BNC  is  “temperature  driven,”  while  in  2DES 
the  transition  is  “density  driven.”  By  this  we  mean  that,  in 
2DES,  the  occurrence  of  metallic  or  insulating  regimes  is 
primarily  determined  by  the  carrier  density  and  it  is  fairly 
independent  of  the  temperature,  so  the  transition  is  defined  by 
a  critical  density.  On  the  other  hand,  in  h-BNC,  both  metallic 
and  insulating  regimes  can  be  observed  for  all  carrier  densities 
studied  in  our  work,  and  the  transition  is  characterized  by  a 
critical  temperature  separating  the  two  regimes. 
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A  second  important  difference  is  that,  contrary  to  2DES,  in 
h-BNC  doping  with  either  electrons  or  holes  seems  to  suppress 
the  metallic  state.  In  other  words,  the  range  of  temperatures  in 
which  the  system  is  metallic  decreases  with  doping.  This  is  a 
novel  and  anomalous  behavior. 

The  robust  nonmonotonic  behavior  of  R(T)  for  all  doping 
levels  studied  in  h-BNC  led  us  to  picture  this  system  as 
described  by  the  coexistence  of  metallic  and  insulating 
conducting  channels,  as  described  in  models  for  the  IMT  such 
as  in  manganites. 40-42  This  picture  may  also  be  adequate  for  a 
composite  system  like  h-BNC.  If  we  assume  these  mechanisms 
act  independently,  we  can  formulate  a  phenomenological 
parallel-resistor  model  to  describe  the  R(T)  curves.  As 
depicted  schematically  in  Fig.  16, 

R(T)  «  — - —  +  - 

L  Rm(T)  Ri(T) 

where  Rm(T )  and  Ri(T)  are  the  metallic  >  0)  and  insu¬ 
lating  (^  <  0)  components  of  the  resistance,  respectively. 

In  the  metallic  regime,  T  rpeak,  then  R(T)  ~  Rm(T).  We 
find  that  Rm{T)  has  the  empirical  exponential  form  also  typical 
of  2DES43 

Rm{T)  =  R0  +  Ri  exp (~Tm/T),  (2) 

for  which  no  final  theory  exists  as  yet,5  with  different 
proposals  as  diverse  as  temperature-dependent  screening,44 
percolation,45  charging  of  traps,46  and  an  analog  of  the 
Pomeranchuk  effect.47 

In  the  insulating  regime,  T  »  Tpeak,  R(T )  -  Ri(T),  the 
conductance  in  2DES  that  is  often  associated48,49  with  either 
Mott50  or  Efros-Shklovskii32,33  (ES)  variable-range  hopping 
(VRH)  between  localized  states.  Our  data  are  better  described 
by  the  ES-VRH  model,  also  known  as  VRH  under  a  Coulomb 
gap: 

Ri(T)  =  R2exp[(Ti/T?'2].  (3) 

Figure  16  shows  a  perfect  agreement  between  this  parallel- 
resistor  model  and  the  experimental  data.  The  best-fit  param¬ 
eters  for  the  90  at.  %  carbon  sample  are  Ro  =  0.2  MQ,  Ri  = 
87  MQ,  R2  =  0.1  M£2,  Tm  =  2.23rpeak,  and  7}  =  2irpeak- 

The  doping  (or  gate  potential)  dependence  of  resistance  is 
automatically  built-in  through  the  proportionality  between  7} , 
Tm,  and  Tpeak,  which  is  a  direct  consequence  of  the  single¬ 
parameter  scaling  of  the  data.  The  low-temperature  limit  of 
the  resistance  is  independent  of  doping,  with  R(T)  ~  Rq  ~ 
0.2  MQ.  This  we  attribute  to  scattering  off  the  insulating  BN 
nanodomains,  since  it  agrees  qualitatively  with  the  resistances 
obtained  from  NEGF  calculations  of  the  transmission  function 
through  several  model  h-BNC  systems,  as  described  in  Sec.  III. 

Our  data  analysis  allows  for  a  clear  identification  of 
the  analytical  expressions  associated  with  metallic-  and 
insulating-type  electronic  conduction.  We  start  by  analyzing 
the  insulating  behavior.  Figure  17  shows  the  R(T)  data  for 
—40  V  gate  voltage  in  the  insulating  region  (T  >  Tpea k)  at 
zero-magnetic  field.  In  the  figure,  ln(R/Rpeak)  is  plotted  against 
(r/rpeak)_x,  so  if  the  data  follow  a  straight  line  it  is  an 
indication  that  the  empirical  law  R(T)  ~  exp[(T/7^)_x]  is 
obeyed  (7)  being  proportional  to  Tpeak). 


FIG.  17.  (Color  online)  Data  fitting  of  the  insulating  portion  of 
the  R(T)  data.  Black,  red,  green,  and  blue  curves  correspond  to  x 
values  of  1, 1/2, 1/3,  and  1  /4,  respectively.  The  ES-VRH  mechanism 
(x  =  1  /2)  gives  the  best  fit  of  our  data  to  a  straight  line. 

Different  values  of  x  correspond  to  analytical  expressions 
commonly  associated  to  different  physical  mechanisms  for 
conduction  in  insulating  materials.  For  instance,  x  =  1  cor¬ 
responds  to  Arrhenius  activation,  which  would  imply  the 
existence  of  a  gap;  x  =  1/2  is  the  Efros-Shklovskii  variable- 
range  hopping  (for  any  dimensionality);  and  jc  =  1/3  and 
x  =  1  /4  correspond  to  Mott  variable-range  hopping  for  two 
and  three  dimensions,  respectively.  As  indicated  by  the  straight 
lines  (guides  to  the  eye),  the  Efros-Shklovskii  hopping  is 
clearly  favored  in  a  larger  temperature  range  with  respect  to 
the  other  fitting  forms. 

Efros-Shklovskii  variable-range  hopping  (ES-VRH),  also 
known  as  VRH  under  a  Coulomb  gap,  is  a  well-known 
conduction  mechanism  in  insulators32,33  and  it  requires  the 
following  ingredients:  (1)  a  distribution  of  localized  states 
around  the  Fermi  level  and  (2)  long-range  Coulomb  repulsion 
between  electrons  in  these  localized  states.  These  conditions 
will  lead  to  a  Coulomb  gap  in  the  electronic  DOS  that  in  two 
dimensions  is  actually  a  “soft  gap,”  that  is,  a  linear  dip  in 
the  DOS  precisely  at  the  Fermi  level.  VRH  between  states  at 
such  nonconstant  DOS  leads  to  the  ES-VRH  expression  for 
the  resistivity  R(T)  ~  exp[(7y  7})1/2],  with 

e 2 

Ti=P  — ,  (4) 

KCl 

where  is  a  constant,  a  is  the  localization  radius  of  the 
states,  and  k  is  an  effective  dielectric  constant.  The  doping 
dependence  of  7},  that  decreases  for  both  electron  or  hole 
doping,  is  consistent  with  a  larger  screening  of  the  Coulomb 
interactions  upon  enhancing  the  number  of  carriers  in  the 
percolating  graphene  network  around  the  BN  islands. 

We  now  turn  our  analysis  to  the  metallic  behavior.  Figure  1 8 
shows  the  R(T)  data  for  +28  V  gate  voltage  in  the  metallic 
region  (T  <  Tpea k)  at  zero-magnetic  field.  In  the  figure, 
In [(R  —  RoV^peak]  is  plotted  against  Tpeak/  T.  The  data  follow 
very  closely  a  straight  line,  indicating  that  the  empirical 
law  R(T)  =  Ro  +  Ri  exp (—Tm/T)  is  obeyed  (Tm  also  being 
proportional  to  Tpea k). 

The  model  described  above  implies  the  coexistence  of 
percolating  graphenelike  electronic  states  that  are  responsible 


(1) 
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FIG.  18.  (Color  online)  Data  fitting  of  the  metallic  portion  of  the 
R(T)  data. 

for  the  metallic  channel  and  a  distribution  of  localized 
electronic  states  through  which  VRH  occurs.  These  states  must 
coexist  at  the  Fermi  level  for  all  considered  doping  levels,  but 
must  also  be  spatially  separated,  so  that  the  two  channels  are 
uncoupled. 

From  DFT  calculations  for  both  N-terminated  zigzag  edged 
nanodomains  and  zigzag  edged  nanoribbons  embedded  in 
graphene,  we  find  indeed  a  coexistence  of  localized  edge 
states  at  the  nanodomain  boundaries  [Fig.  13(a)]  and  extended 
graphenelike  states  [Fig.  13(b)]  near  the  Fermi  level  (separated 
by  about  0.2  and  0  eV,  respectively),  for  the  range  of  doping 
considered  here.  Since  these  embedded  extended  BN  ribbon 
and  small  triangular  BN  domain  calculations  describe  the 
limiting  cases  of  large  and  small  h-BN  nanodomains,  it  is 
reasonable  to  assume  that  both  localized  edge  states  and 
extended  graphenelike  states  are  present  in  the  experimental 
system.  Disorder  and  potential  fluctuations  will  likely  broaden 
the  energy  distribution  of  localized  edge  states,  with  Coulomb 
repulsion  causing  the  insulatorlike  ES-VRH  conduction  as 
opposed  to  metal-like  band  conduction.  Spatial  separation  of 
the  two  channels  will  occur  by  the  simultaneous  percolation 
of  “metallic  regions”  and  “insulating  regions.” 

Finally,  we  comment  on  the  doping  dependence  of  the 
characteristic  insulating  and  metallic  temperatures  7}  and  Tm , 
which  are  both  proportional  to  rpeak-  Since  Ypeak  decreases 
upon  doping  with  either  electrons  or  holes,  so  do  7]  and  Tm . 
This  suggests  that  doping  affects  the  resistance  of  h-BNC  via  a 
common  mechanism  for  both  metallic  and  insulating  regimes. 

As  described  above,  the  decrease  of  7)  with  doping  may  be 
qualitatively  understood  within  the  well-established  theory  of 


ES-VRH,  since  7)  should  decrease  with  increasing  screening, 
which  should  by  itself  increase  with  doping.  However,  the 
decrease  of  Tm  with  doping  does  not  fit  the  phenomenology 
of  the  metallic  phase  of  semiconducting  2DES  nor  any  of 
the  related  models.  In  fact,  it  shows  precisely  the  opposite 
behavior. 

Of  course,  graphene  systems  are  in  many  ways  different 
than  semiconductor  2DES,  not  only  for  their  band  dispersions 
and  DOS  profile  near  the  Fermi  level,  but  also  because  in  doped 
2DES  carriers  are  either  electrons  or  holes,  but  in  h-BNC 
(as  in  graphene),  potential  fluctuations  cause  the  coexistence 
of  electron  and  hole  puddles.25  Therefore,  the  present  results 
represent  a  novel  type  of  IMT  in  2D  electronic  systems,  and 
they  will  likely  open  a  new  pathway  in  the  investigation  of 
such  systems. 
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